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Abstract 

We examine a hot pion gas by including medium modifications of the two-body 
scattering amplitude as well as mean-field effects selfconsistently. In contrast to 
earlier calculations, the in-medium T-matrix is rather close to the free one while 
the mean-field potential agrees well with lowest-order estimates. We also discuss 
the validity of the quasiparticle approximation. It is found that it is reliable for 
temperatures up to ~ 150 MeV. Above this temperature off-shell effects in the pion 
selfenergy become important, especially if the pions are strongly out of chemical 
equilibrium. 
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1 Introduction 



In ultrarelativistic heavy ion collisions several hundred particles are produced in the final 
state. In the midrapidity region at CERN energies of 200 GeV/A most of them are pions 
(with a pion-to-baryon ratio of ~ 6 : 1). The lifetime of the pionic fireball is a few fm/c 
and the freeze-out volume is typically 10 3 fm 3 . The theoretical understanding of such a 
thermal 'pion gas' is currently of great interest in connection with possible signals from the 
quark-gluon plasma (QGP) phase transition 0]. On the other hand a hot and dense Bose 
gas is also interesting from a many-body point of view since one may expect interesting 
correlations associated with the statistics. 

Several aspects of the thermal pion gas have been discussed previously, including mean- 
field effects [[| and medium modifications of the cross section ||. It is the purpose of 
the present paper to give a unified description of these processes by using a reliable model 
for the vacuum nn interaction and requiring selfconstistency at the two-body level. This 
leads to a Brueckner scheme, familiar from the microscopic many-body theory of nuclear 
matter. Throughout the discussion we shall assume that the gas is in thermal but not 
necessarily in chemical equilibrium ({i n ^ 0). The latter seems to be required from fits 
to pion ]?T-spectra, which yield fi^ ~ m n . The assumption of thermal equilibrium, on 
the other hand, seems reasonable from the following simple estimate: at 200 GeV/A the 
freeze-out density n n of pions produced at midrapidity is ~ 0.3/m -3 . With N n ~ 400 this 
yields a freeze-out radius Rf ~ 7 fm. Taking an average nn cross section (a) = 15 mb gives 
a mean free path A = l/n 7T (a) ~ 2.2 fm and hence the mean number of collisions is ~ 3. 
Thermalization is further corroborated in a scenario where initially a quark-gluon plasma 
(QGP) is formed. Here QCD string breaking models as well as parton cascade simulations 
also yield thermal equilibration which should survive during the hadronization. 
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2 The Vacuum tttt Interaction 



The starting point for our description of the hot pion gas is the vacuum tttt interaction 
model of ref. || which is based on meson exchange. Here the basic meson-meson inter- 
action is constructed from an effective meson Lagrangian with phenomenological form 
factors at the vertices. From these vertices two-body pseudopotentials are constructed in- 
cluding the most important s- and ^-channel meson exchanges. Among these the ^-channel 
p exchange between tttt states and the s-channel p pole term will be most important for our 
discussion. We employ the Blankenbecler-Sugar (BbS) reduction || of the 4-dimensional 
Bethe-Salpeter equation maintaining covariance 0. The irir T-matrix for given angular 
momentum J and isospin I is obtained as 



+ 



/ dk e Aul V^(Z, Ql , k) G^(Z, k) T^(Z, k, q 2 ) , (1) 
Jo 



where k = \k\ etc. ; Z is the CMS energy and G® n (Z, k) the vacuum two-pion propagator 
in the CMS frame with pions of momenta k and —k (the Lohse et al. model also contains 
coupling to the KK channel which has been omitted for brevity in eq. (|XJ) but which is 
included in the calculation). In the BbS form the two-pion propagator is given by 

with ojk 2 = k 2 + ml. This model gives a good description of the phase shifts and inelas- 
ticities up to ~ 1.5 GeV which is more than sufficient for our purposes. 



3 Selfconsistency 

The most obvious medium modification of the tttt scattering in the gas surrounding is 
a change in momentum weight of the intermediate two-pion propagator, first studied in 
ref. 0. Identifying the CMS frame with the thermal reference frame one has 

r (7 u t\- 1 ( 1 + Vk(^,T)) 

Lr nw {Z,, K, p w , 1 ) — — — "2 , [O) 
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where fk = (exp[(ujk — [i w )/T] — is the thermal Bose factor and fi n the chemical 
potential. The identification of the CMS frame with the thermal frame simplifies our 
calculations considerably (allowing for a relative velocity between the two frames we find 
that it effectively acts like a change in fi n which is a parameter for us anyway). At fixed 
T the pion chemical potential fixes the density via 

n w (fi w ,T) = g w J fqin^T) , (4) 

where g w = 3 is the isospin degeneracy factor. We use a temperature range of 100-200 
MeV. A temperature of 100 MeV is roughly the lowest temperature from thermal fits 
to pr-spectra at the AGS. On the other hand 200 MeV should be an upper limit for 
purely hadronic models since one expects the phase transition into a QGP around this 
value. While thermal equilibrium seems a reasonable assumption it is not clear whether 
chemical equilibration is reached during the evolution of the pion gas. Indeed, fits of the 
CERN pT-spectra could be improved with ^ ~ 125 MeV J4|. Hence we discuss both 
cases fi n = and fi n = 125 MeV. At the same temperature the latter gives a higher pion 
density reaching a maximal density n n ~ 0.7 fm~ 3 at 200 MeV. 

There is a second effect which needs to be considered [0 . The tttt interaction introduces 
a pion selfenergy T, n ('mean field') which changes the single-pion dispersion relation 

u]\ = m 2 n + k 2 + E n (uj k , k;ii n ,T) (5) 

as a function of density and temperature. In terms of the forward-scattering amplitude 
M n7T , is expressed as [Q, [| : 

E^A^T) = J tJ^^ f P (^,T) MUk {A) J A) ) ■ (6) 
Relating the forward-scattering amplitude to the T-Matrix as 

M^{E cms ) = (2tt) 3 E 2 cms T nn (E cms ) , (7) 
where E cms = y /s=[(p+ ky (p + k)^} 1 / 2 is the CMS energy of the two colliding pions, 
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one can transform the selfenergy expression into 

OO Emax 

E n (u,k;fJ,n,T) = ^ J dp -j- f p (fjLn,T) J dE cms E^ ms T^{.E cms ) . (8) 

P F, - 

Here we have restricted ourselves to the on-shell T-Matrix neglecting the dependence on 
the total momentum P = k + p of the pair. Thus the energy integration bounds are given 

as 

Emax/min = + + 2uU p - k 2 - p 2 ± 2kp) 1/2 . (9) 

For the forward-scattering T-Matrix we take the spin-isospin weighted sum including 
partial waves up to J = 2: 

TUEcms) = ±- E (2 J + 1) T^(£ cms ) , (10) 

which saturates the cross section in the relevant energy range. 

It should be noted that the rate = —21-mL^I (2(k 2 + m 2 ) 1 ^ 2 ) deduced from eq. @ is 
not entirely consistent with that obtained from the collision term of the bosonic Boltzmann 
equation. A correct account of the Bose correlations of the two interacting pions will 
modify the occupancy factor in (|8]) ||. It turns out, however, that this is a small effect 

Hi- 
lt is now evident that the pion selfenergy and the T-Matrix should be combined in a 
selfconsistent Brueckner scheme as indicated in Fig. 1. This implies that the selfenergy 
(8) is to be calculated from the in-medium T-Matrix which, on the other hand, should be 
obtained from the in-medium 27T propagator including the pion selfenergy: 

G 7r7T (Z,k;^,T) = (l + 2f k ( f i n ,T)) [ ^ D v (u, k) D V {Z - u, k) , (11) 

J 27T 

where 

D„(u, k) = [to 2 -ml-k 2 - E w (w, k; ^ T)]~\ (12) 
Together with the in-medium scattering equation, 
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equations (8), (11) and (12) define the selfconsistency problem (see also Fig. 1). It should 
be noticed that the pseudopotential V nn remains unchanged. 

In the following we will discuss two different methods of calculating G^. 



4 Quasiparticle Approximation 

The quasiparticle approximation (QPA) is valid if the lifetime is long or more precisely if 
the quasiparticle energy 

e k = {ml + k 2 + Re^{e k ,k)fl 2 (14) 

is much larger than its decay width IV In this case the energy-dependence of ReT,^ is 
expanded to first order around the 'quasiparticle pole' e k as 

ReE n (u, k) « ReE n (e k , k) + 9 ^j?' k) \e k - 4) ■ ( 15 ) 

One can then perform the folding integral (11) analytically which gives 

G^{Z, k; Hn, T) = - — — — (16) 



with 



dReH^Uj k) . , 1 
IhJ 2 

Col = e t + * z k M T (et, k) quasipion dispersion relation . (17) 



z k = (1 — le fc ) 1 the polestrength 



Together with this approximation eqs. (8), (11) and (13) are solved by iteration starting 
from the free pion dispersion relation and the vacuum T-Matrix. We keep the pion density 
n n fixed during the iteration by readjusting in each step (the final [i w differs from the 
starting value /i^ -* by a small amount). The selfconsistent results are shown in Fig. 2. 
Already for chemical equilibrium (//^ = 0) we find a considerable reduction of the peak 
values in ImT JI as compared to the vacuum case in both s- and p-wave. This is in 
agreement with the lowest-order results from refs. |Tl| . The reduction is mainly caused 
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by the Bose factors (1 + 2f) leading to a stronger weighting of the lower pion energies. 
For the same reason the near threshold region shows an enhancement over the vacuum 
T-Matrix especially in the s-wave. We define the pion 'optical potential' as 

which is shown in the lower part of Fig. 2. In agreement with refs. @, [8| we find attraction 
in ReV n for low momenta. The selfconsistent potential, however, is significantly weaker 
than in lowest order. The maximum in ImV n is due to formation of the p-resonance, 
whereas the non-zero values at k = arise from s-wave interaction with thermally moving 
pions. 



5 Off-shell Integration of G nn 

To check the reliability of the QPA we perform the same calculations as described in 
the previous section, but the quasiparticle two-pion propagator is now replaced by a 
numerical integration of (11) accounting for the full off-shell properties of S 7r (a;, k). Using 
the symmetry relation 

E„(-cjJ, k) = E 7r (o;, k) , (19) 
the in-medium propagator (11) can be written as 

G^(Z,k;^,T) = (l + 2/ fc (/i 7r ,T)) / — D v (u,k) D V [\Z - u\,k) . (20) 

J 7T 

Z/2 

Fig. 3 shows the selfconsistent results for a chemically equilibrated pion gas (/j,^ = 0). At 
T = 125 MeV the results coincide with the QPA within a few percent. For T^150 MeV the 
deviations become larger: the T-Matrices in s- and p-wave are now enhanced compared 
to the vacuum curve for most of the energy range (e.g. the peak value of the p-resonance 
increases by ~ 30% at T = 200 MeV). This clearly must be an off-shell effect. From 
the lower part of Fig. 3 one can conclude that a first-order expansion of ReY^^ around 
efc does not describe the energy dependence correctly; in addition Y k jt k ~ 0.6 around 
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k = 300MeV/c, i.e. the quasiparticle lifetime becomes very short and renders the QPA 
invalid. The potentials V n (k) are now very close to the lowest-order results. As compared 
to the QPA a considerable amount of attraction is restored in the low-momentum region 
of ReV- K {k). We also investigate a scenario with finite chemical potential fi^ = 12hMeV 
as suggested by thermal fits to the SPS ]>r-spectra ||. The s-wave ttti interaction now 
shows a strong accumulation of strength near threshold (Fig. 4) which might lead to quasi- 
bound 7r + 7r~ pairs as suggested in ref. [n|. The latter possibility deserves further study. 



In the p-wave the changes of the p-resonance are appreciable: at highest temperatures the 
width decreases considerably which might be detectable via dilepton pairs coming from 
the midrapidity region. 

6 Summary 

Based on a selfconsistent Brueckner theory we have presented a numerical analysis of a hot, 
interacting pion gas in thermal equilibrium, which may be realized in future experiments 
at RHIC or LHC (at ^/i > lOOGeV/A). The selfconsistent Brueckner scheme accounts 
for statistical (Bose factors) as well as dynamical (selfenergy) modifications of the pion 
propagation in the gas. Using a realistic model for the vacuum interaction we have 
solved the non-linear problem by iteration. We have compared results from the QPA 
to those taking full account of the off-shell behavior in the pion selfenergy. It is found 
that the QPA breaks down for temperatures T > 150MeV and finite chemical potentials, 
i.e. high pion densities. In the full calculation we find a slight increase of the in-medium 
T- matrix - in contrast to refs. 0, where only statistical modifications have been taken 
into account. The main cause of this enhancement is attributed to strong off-shell effects 
in the pion selfenergy. The single-particle potentials show considerable attraction for low 
pion momenta, although quantitatively somewhat less than our lowest-order results. As 
was shown recently by Koch and Bertsch jl2| , such an attraction is not able to explain the 



low-pr enhancement in the pion spectra of current experiments. However, as suggested 
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by numerical simulations of the bosonic transport equation [13|, an increase of the in- 
medium T-Matrix, as found in our calculations, is likely to ensure thermalization of the 
pionic fireball and hence to produce an excess of low-j>r pions through the (1 + /) factors 
in the collision integral. This issue and the impact of finite baryon density in the pion 
gas will be addressed to a future publication. 
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FIGURE CAPTIONS 



Fig. 1 Set of selfconsistent equations for an interacting pion gas: 
upper part: Dyson equation for in- medium pion propagation 
lower part: in- medium tin T-Matrix equation. 

Fig. 2 Selfconsistent pion gas at ^ = in QPA: 

upper part: imaginary part of the rnr T-Matrix in the JI=00- and JI=ll-channels for 
several temperatures (long-dashed lines: T = 125 MeV, short-dashed lines: T = 150 
MeV, dotted lines: T = 200 MeV) 

lower part: corresponding pion potentials (short-dashed lines: T = 150 MeV, dotted 
lines: T = 200 MeV; the dashed-dotted and dashed-double-dotted lines show the lowest- 
order results calculated with the vacuum T-Matrix for T = 150 MeV and T = 200 MeV, 
respectively). 

Fig. 3 Selfconsistent pion gas at fi^ = with full off-shell integration of G nn (eq. (20)): 
upper part: imaginary part of the nn T-Matrix (see Fig. 2) 

lower part: real part of the pion potentials (see Fig. 2) and real part of the off-shell pion 
selfenergy at T = 200 MeV. 

Fig. 4 Imaginary part of the selfconsistent T-Matrix with full off-shell integration of 
G wn for = 125 MeV (long-dashed lines: T = 100 MeV, short-dashed lines: T = 150 
MeV, dotted lines: T = 175 MeV, full lines: free space). 
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